Outer-gap vs. Slot-gap Models for Pulsar High Energy Emissions: 

The Case of the Crab Pulsar 



Kouichi Hirotani 

ASIA A/National Tsing Hua University - TIARA, 
PO Box 23-141, Taipei, Taiwan^ 
hirotaiiiOtiara. sinica. edu.tw 



ABSTRACT 

We analytically examine the capabilities of rotation-powered pulsars as the sources of gamma- 
rays and show that their phase-averaged gamma-ray flux is proportional to the product of the 
spin-down flux and the gap trans-field thickness cubed irrespective of the emission models. Ap- 
plying the scheme to the Crab pulsar, we demonstrate that the outer-gap model reproduces the 
observed GeV fluxes and that the slot-gap model reproduces at most twenty per cent of the 
observed fluxes because of the small trans-field thickness. An implication on the relationship 
between the gamma-ray and the spin-down fluxes is discussed. 

Subject headings: gamma-rays: observations - gamma-rays: theory - magnetic fields - methods: analyt- 
ical - pulsars: individual(Crab) 



1. Introduction 

The launch of Fermi Gamma-ray Space Tele- 
scope will soon open a new era for the studies of 
rotation-powered pulsars. The unprecedented sen- 
sitivity and spectral resolution of the Large Area 
Telescope (LAT) aboard Fermi Space Telescope 
will allow for detailed studies of particle accelera- 
tion and radiation in rotating neutron-star (NS) 
magnetospheres. To make the best use of the 
power of LAT observations, we need the most so- 
phisticated model with minimum assumptions. 

In all the pulsar emission models, e~'s and/or 
e"*''s are accelerated by the magnetic-field-aligned 
electric field, E\\, to radiate photons in the open 
zone (fig. 1 in Hirotani 2008, hereafter H08) 
mainly via synchro-curvature process. In polar- 
cap (PC) models, emission takes place within 
several NS radii above a PC surface (Arons & 
Scharlemann 1979; Daugherty & Harding 1982, 
1996). However, such a low-altitude emission pre- 
dicts too small beam size to produce the observed 
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wide pulse profiles. Therefore, extending the orig- 
inal idea by Arons (1983), Muslimov and Harding 
(2003, hereafter MH03; 2004a, b) and Dyks et 
al. (2004, hereafter DHR04) sought the possibil- 
ity of a wide hollow cone of high-energy radiation 
due to the flaring of field lines. They proposed 
the slot-gap (SG) model, in which emission takes 
place along the last-open field lines. Recently, 
Harding et al. (2008, hereafter HSDF08) demon- 
strated that the SG model reasonably reproduces 
the Crab pulsar phase-resolved spectrum (Fierro 
et al. 1998; Kuiper et al. 2001; Nolan et al. 1993). 

The outer-gap (OG) model gives an alternative 
possibility (Cheng et al. 1986a, b; Romani 1996; 
Cheng & Zhang 1996; Hirotani 2006a, b, hereafter 
H06a, b). It differs from the SG model in the 
following ways: • An OG extends between the 
null surface on which the magnetic field becomes 
perpendicular to the rotation axis and the light 
cylinder on which the plasma co-rotational veloc- 
ity would exceed the speed of light, c. A SG ex- 
tends between the PC surface and the light cylin- 
der. • For Crab, an OG occupies more than 10% 
of the open magnetic fluxes with a large electro- 
static potential drop, A^E" ~ 10^^ V, whereas a 
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SG occupies several % of the open fluxes with 
Avf ~ 10^"^ V. • Pair production copiously takes 
place in OGs, whereas it is negligible in SGs. 

The purpose here is to explore further into the 
two models. We examine their common emission 
properties in § [2l and separately consider the OG 
and SG models in §§[3]and[31 §[5]is for discussion. 

2. Gamma- ray Flux 

In this letter, we concentrate on the phase- 
averaged spectrum of pulsar magnetospheric emis- 
sions, sacrificing the examinations of light curves 
and phase-resolved spectra. In this context, we 
can neglect the aberration of photon propagation 
directions and the time-of-flight delay due to dif- 
ferent emission altitudes from the NS, because 
these two relativistic effects do not change the to- 
tal number of photons to be detected. Although 
this thought experiment does not describe any re- 
alistic pulsar emissions, it significantly reduces the 
calculation of photon propagations and gives the 
correct phase- averaged spectrum. 

To investigate the upper hmit of photon fluxes, 
we neglect photon absorptions. Then the radia- 
tive transfer equation gives the specific intensity 
Jj/ « 2bgcju, where 6 « 7"^ denotes the emission 
beaming angle, 7 the Lorentz factor of e^'s, Qc 
the local curvature radius of the magnetic field 
line, jj/ the emission coefficient, and 2bgc the dis- 
tance interval from which the photons are de- 
tected by the observer (fig. 6.2 in Rybicki & Light- 
man 1979). Giving the emission coefficient as 
> w N{dP/diy)/{TTb'^), we obtain 

TT dv 

where N denotes the spatial density of e~'s or 
e"'"'s, and dP/dv the radiation power per particle. 

At each magnetic azimuth Lp^ on the PC sur- 
face, we parameterize the field-line footpoint with 
their magnetic colatitude 0* measured from the 
magnetic axis. We define that the primary pho- 
tons are emitted only along the field lines thread- 
ing the PC surface with ef'"" < 9^ < Of'"'', where 
the upper boundary and the last-open field lines 
are designated with 0™™ and 9f^^, respectively. 
We further introduce the dimensionless gap trans- 
field thickness, /i,„ = {Of'^'' ~ 9T'')/dT''. In this 
letter, we assume that is constant for both Lp^ 



and s (distance along the field line) for simplicity. 

Let us introduce the gap meridional thickness 
A2 that represents the distance between the last- 
open field line and the upper boundary measured 
perpendicularly to the field line. Then, the ob- 
server detects emissions from the magnetic-flux 
cross section /S.A k, Az x 2forsin6', where rsn\9 
refers to the distance from the magnetic axis, 
and 2b the azimuthal full opening angle of the 
points from which the photons propagate towards 
the observer. For an aligned rotator, the az- 
imuthal length from which the photons propagate 
towards the observer is given by 26rsin0; thus, 
/S.A = Az X 2br sin 9 holds. For an oblique rotator, 
in the outer magnetosphere, toroidal expansion of 
the field line flux tubes is similar to an aligned 
case; thus, AA « Az x 26rsin6' approximately 
holds. Since Az « 2h^vji^c[r /vji^q)"^ csc9{l + 
3cos^ 0)^^/^ holds for a dipole fleld, we obtain 

AA » 2bK,vjlc{B,/B){r,/vDi^c)\ (2) 

where r denotes the distance from the NS center, 
9 the magnetic colatitude, -culc = c/fi the light 
cylinder radius, f2 the NS angular frequency. The 
magnetic field strength B is evaluated at the emis- 
sion point and related with its surface value -B, by 
= (n/r)Vl + 3cos2 61/2. 
The photon energy flux can be computed by 

vF, = uI,{^A/d^), (3) 

where d is the distance to the pulsar. Substituting 
equations H]) and ^ into equation ([3]), we obtain 

vFi, « -h^- ngcV— I — , (4) 

TT 27rce dv \ d / t^lc 

where e is the charge on the positron; the dimen- 
sionless particle density per magnetic flux tube, 
n = (2-Kce/VlB)N , becomes approximately cosai, 
where ai is the magnetic inclination with respect 
to the rotation axis. The quantity ^lB^/{2TTce) 
expresses the typical Goldreich-Julian (GJ) parti- 
cle number density at the PC surface (Goldreich 

6 Julian 1969). It follows from equation ^ that 
the photon flux vFi, does not depend on b. 

For saturated e~'s or e+'s, electrostatic force 
balance, ei?|| = 2e'^j'^/{3g^), gives the terminal 
Lorentz factor 

/ ^ 1/2 / p s 1/4 

7 = 2.20xlO'5l]-^/M-^) (-^) , 

(5) 
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where ^2 = f7/(102 rads"^) and = -(B/B) ■ 
V\l/. The non-corotational potential "if is given by 
the inhomogeneous part of the Maxwell equations, 

- V^* = 47r(p - pGj) = i2nB,/c)K, (6) 

where p and pcj denote the real and GJ charge 
densities, respectively, and Bz the magnetic field 
component projected along the rotation axis. The 
effective charge density, pcj — p, is parameterized 
by K, which is a function of position. In the New- 
tonian limit, we would have k = 1 for a vacuum 
gap while k < 1 for a non- vacuum gap. 

Since e^'s are ultra- relativistic in the gap, 
the primary emission is dominated by the pure- 
curvature component. Thus, 



,— 

dv 



3a/3 06^74 „ 
47r p2 



K^/mdi (7) 



peaks at X = v/vc — 1.318 with the maximum 
value x'^ K5/3{£_)d^ — 0.6826, where is 
the modified Bessel function of 5/3 order, and 
Vc = (3/47r)(c7'^/pc) the characteristic frequency 
of the emission. The electrostatic force balance, 
(2/3)eV/e? = e£^|h gives 



dP 9V3 ^ ^ 
— — — — — cehnx' 



dv 



87r 



K^/^dt (8) 



For a thin gap (/i,„ ^ 1), equation (68) in H06a 
gives 



hi, ( 



4 V^LC 



b 9{-.BJB) 



where s is the distance along the field line. For the 
field lines curving away (or toward) the rotation 
axis, 1.0 < -d{BjB)/dls/wi^c) < 1-5 (or 1.5 < 
—d{Bz/B)/d{s/vui^c) < 2.2) is the typical range. 

Substituting B» = 2/^/rJ and i^5/3(C)'^^ = 

0.68, and combining equations ([U, ^ and we 
finally obtain the peak flux. 



peak 

w 0.0450//i^K 



2^(4 1 



d2 



(10) 



where p, denotes the magnetic dipole moment and 
I g, 1 d{-KBjB) 



f 



0.7 0.5ci7Lc 1-5k 9(s/n7Lc) 



(11) 



is close to unity (see also Zhang & Cheng 2003). 
If we apply ^5/3(0^^ = 1.61, we 



also obtain the integrated flux, F^dv = 
2.36(i/i^i,)pcak- Note that the factor p'^Q.^ / {c^d^) 
is proportional to the spin-down flux (§[5]). 

It should be emphasized that equation PH]) 
estimates the upper limit of the ]ihase- averaged 
flux. For example, for an axisymmetric gap in an 
aligned rotator (ai = 0°), the phase-averaged flux 
equals equation (jlOp . which is constant during the 
whole NS rotation. Since equation (|10p holds par- 
ticularly well for an aligned rotator, and since the 
phase-averaged flux decreases with increasing a\ 
(e.g., MH03 and DHR04), equation ^ gives the 
upper limit for general ai. The actual flux will 
be a few times less than equation (ITUl) . because 
only the field lines in a limited azimuthal range 
are active, and because the pcj at the PC surface 
decreases with increasing aj. In § 31 we confirm 
this by comparing with numerical results. 

The difference between the OG and the SG 
models comes into equation (jlOp through hm_, 
and the assumed p. In §§ [3] and IH we apply 
this equation to these two models, considering the 
brightest spin-down-flux pulsar, the Crab pulsar, 
assuming d = 2 kpc. 

3. Outer-gap Model 

First, let us apply equation pO)) to the OG 
models. Using the vacuum (k = 1.0) OG mod- 
els of Cheng et al. (2000), Takata et al. (2008) 
and Tang et al. (2008), which proposed h^a ~ 0.11 
(i.e., / « 0.11 in their notation), we obtain 



(J^i^^.)peak « 6.58 X lO^V^; 



0.11/ V3.8/ 



with E\\ « 2.55 x 10^ V m 

-130, 



MeVs~^ cm"^ (12) 



, 7 « 2.0 X 10^, where 
pzo = ^/(lO^Gcm^). The vF^ flux peaks at 
1.3/iz/c ~ 4.1 GcV. These results are consistent 
with their vacuum OG models and with the ob- 
served phase-averaged spectrum, where the pho- 
ton flux around GeV is dominated by the pri- 
mary curvature component rather than the re- 
processed synchrotron-self-Compton one (Takata 
& Chang 2007). Exactly speaking, the vacuum 
OG model is not self-consistent electrodynami- 
cally, because it assumes a vanishing charge den- 
sity while it adopts the GJ flux of radiating e+'s. If 
one instead adopts the actual e+ flux in a (nearly) 
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vacuum OG, the predicted 7-ray flux would be a 
few orders of magnitude less than equation (|12p 
(left panel of fig. 8 in H06a). That is, equa- 
tion p^ . and hence the traditional OG models 
reproduce observations only phenomenologically. 

Secondly, let us consider the non-vacuum OG 
solution, which consistently takes the pair pro- 
duction in the pulsar magnetosphere into account. 
Using K ~ 0.3 and /im « 0.14, we obtain 

MeVs^^cm"^ (13) 

which is consistent with the numerical results for 
Crab (figs. 8 & 9 in H08). In short, the OG model 
reproduces the observed flux around GeV. 

4. Slot-gap Model 

Thirdly, let us apply equation (fTO]) to the SG 
model, which adopts h^a < 0.05 (i.e., A^sG < 0.05 
in HSDFOS's notation) and k « -(0.1 + 0.2s/roLc) 
for the Crab pulsar. At the outer part (s « 
O.Gtulc), where emissivity is large, we obtain 

(^F.)pcak « 6.02 X 10- V— ( —V (^Y 
^ 0.22 V0.05 J \ 8 J 

MeVs-^cm-^ (14) 

with l^ii I w 2.4 X lO'^ V m-^ and 7 « 1.1 x lO'^. 

If we adopt the same parameter set = 0.04 
and — 11 as HSDF08, and if we adopt \E\^ \ « 
7 X lO^Vm"^, which is derived from the poten- 
tial drop 0.5 X 1.3 x 10^^ V in the higher alti- 
tudes SG (HSDF08), equations (H and JS]) give 
(j^J^i.)peak ~ 1-9 X 10-5MeVs"^cm-2. This an- 
alytical prediction can be conflrmed by a numer- 
ical computation of the 3-D SG model (H08). If 
we adopt the same parameter set as HSDF08 and 
if we adopt the same as their equation (4) 
for ai = 45°, we obtain the photon map as flg- 
ure [1] which shows a caustic emission from the 
higher altitudes by virtue of the exclusion of the 
strong lower-altitude emission. Specifying the ob- 
server's viewing angle with respect to the rota- 
tion axis, and integrating over the entire NS ro- 
tation, we obtain the phase-averaged spectrum 
(fig- ©J which lies much below the observed value. 
Since \E\i\ w 7 x lO^Vm"^ gives 7 w 8.3 x 10^ 



(eq. [5]), we interpret that HSDF08, who adopted 
7 ~ 2 X 10^, overestimated (j^-F'i/)pcak by the factor 
(2/0.83)'' - 33. 

Equation (|14p could give the observed flux if 
^J'30 > 17. However, such a large ^30 is not al- 
lowed in recent analyzes of force-free electrody- 
namics. For example, Spitkovsky (2006) derived 
E = inilKi{fi'^rt'^/c^){l + sm'^ ai), where / denotes 
the NS moment of inertia and fl the temporal 
derivative of il. Imposing £45 = //(lO^^ gcm'^) < 5 
(with = 16 km and M < 2.5Mq; e.g., Lattimer & 
Prakash 2000), we obtain ^30 < 3.10\/?4^< 7.07 for 
Crab. Thus, even HSDFOS's value, ^30 = 11, may 
be a little too large. In short, equation gives 
the conservative upper limit and the SG model can 
explain at most 20% of the observed flux. 

5. Discussion 

If we assume that the spin down follows the 
dipole radiation formula, E = 2iJ?il'^ siii^ ai/{3c^), 
equation pO|) becomes 

E 707 

(j^F.)pcak« 0.0954//i^K— — (15) 
d'' sm ai 

where E/(P denotes the spin-down flux at the Sun. 
It was, therefore, natural that the largest spin- 
down-flux pulsars were preferentially detected 
with the Energetic Gamma Ray Experiment Tele- 
scope (EGRET). The same tendency can be pre- 
dicted for Fermi Space Telescope. 

For young pulsars, both /i„i and k are less than 
unity (e.g., hm ~ 0.1 and k « 0.3 for Crab) owing 
to the copious pair production in the magneto- 
sphere. As a result, F^dv becomes much less 
than E/d"^ as demonstrated by H08. For middle- 
aged pulsars, on the contrary, hm > 0.5 and k « 1 
holds (Hirotani et al. 2003), leading to an increas- 
ing ratio of F^dv / {E / d"^) with age. 

To explain the observed relationship F^dv oc 
[E/d^Y''^ of pulsed 7-ray emissions (Thompson et 
al. 1994; Nel et al. 1996), it is essential to examine 
the evolution of /im with age. From equation pU)) . 
we can at least state that the index a defined 
by Fi,dv cx {E/d^)°' becomes less than unity, 
because increases with decreasing E, as dis- 
cussed just above. To examine the evolution of 
hm, we must solve the screening of due to the 
discharge of the produced pairs in 3-D pulsar mag- 
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Fig. 1. — 3-D SG prediction: Photon flux per 
phase [MeVs~^ cm^^ deg^^] at distance d — 2 kpc 
from the Crab pulsar as a function of the viewing 
angle and the phase. Only the emission from the 
SG connected to the north pole (i.e., one of the 
two SGs) is depicted for clarity. 
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Fig. 2. — 3-D SG prediction: Crab phase averaged 
spectrum at three discrete viewing angles: 100° 
(sohd), 110° (dashed), and 120° (dotted). Emis- 
sion from both SGs are considered. For observa- 
tional data points, see Kuiper et al. (2001, and 
references therein), and Aliu et al. (2008) for the 
flux at 25 GeV. 



netospheres. In subsequent papers, we shall look 
more carefully into this issue, by simultaneously 
solving equation Q , the Boltzmann equations for 
e^'s, and the radiative transfer equation under 
minimum assumptions. 
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